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LONG  WAVELENGTH  LIMIT  OF  INTERCHANGE  INSTABILITIES 
I.  INTRODUCTION 

A  considerable  eaounC  of  DNA-sponsored  research  has  been  directed  at 
understanding  interchange  instability  processes  in  the  ionosphere.  This  is 
an  iaportant  topic  since  these  instabilities  can  generate  large-scale, 
long-lasting  plasaa  structure  in  the  ionosphere  which  is  deleterious  to 
radar  and  conmunications  systems.  The  immediate  concern  is  to  understand 
structuring  mechanisms  associated  with  a  BANE  in  order  to  obtain  a 
predictive  capability  which  could  aid  the  operations  of  these  systems  in  a 
nuclear  environment.  To  achieve  this  goal,  research  has  also  been  carried 
out  to  understand  ionospheric  structuring  in  a  non-nuclear  environment. 
This  work  involves  the  study  of  natural  structuring  processes  which  occur 
during  equatorial  spread  F  and  in  the  high  latitude  ionosphere,  as  well  as 
structuring  processes  associated  with  artificial  plasma  cloud  releases. 
The  underlying  mechanism  believed  responsible  for  structure  in  all  of  these 
situations  is  an  interchange  instability,  i.e. ,  the  interchange  of  high  and 
low  density  plasma  in  the  presence  of  an  external  force.  The  instabilities 
Investigated  to  date  are  (1)  the  Rayleigh-Taylor  instability  (driven  by  a 
gravitational  force),  (2)  the  E  x  B  gradient  drift  instability  (driven  by 
an  ambient  electric  field  or  neutral  wind),  and  (3)  the  current  convective 
instability  (driven  by  a.  current  parallel  to  the  magnetic  field).  The 
extensive  research  of  these  instabilities  has  entailed  both  theoretical  and 
computational  analyses.  The  theoretical  studies  have  focussed  on  the 
fundamental  linear  and  nonlinear  behavior  of  the  instabilities,  while  the 
computational  studies  have  primarily  dealt  with  the  gross  evolution  of  the 
plasaa  in  the  nonlinear  regime.  For  reviews  of  this  work  see  Ossakow 

(1979),  Fejer  and  Kelley  (1980),  and  Ossakow  et  al.  (1982). 

Ifanmerlpt  spptowd  Nov mabar  1,  1983. 
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A  key  feature  of  the  scmccurlag  of  ionospheric  plasmas  la  both 
nuclear  and  non-nuclear  environments  is  that  large-scale  structure 
develops.  It  has  been  observed  experimentally  and  in  computer  simulations 
that  the  wavelengths  of  the  irregularities  that  develop  (i.e.,  strlations) 
are  generally  comparable  to  the  scale  length  of  the  density  gradient  in  the 
plasma.  However,  the  majority  of  theoretical  analyses  performed  thus  far 
have  considered  irregularity  wavelengths  mnch  smaller  chan  the  density 
gradient  scale  length,  i.e.,  local  theory  has  been  assumed  which  implies  IcL 
»  1  where  k  is  the  wavenumber  and  L  is  the  scale  length  of  the  density 
gradient.  Recently,  the  opposite  limit  has  been  considered  by  Huba  and 
Zalesak  (1983)  in  a  study  of  the  E  x  B  instability  and  by  Huba  (1983)  in  a 
study  of  the  current  convective  instability.  The  crux  of  these  studies  is 
that  a  nonlocal  mode  equation  is  derived  and  solved  analytically  for  each 
instability  in  the  long  wavelength  limit,  i.e. ,  kL  «  1.  Thus,  they 
provide  insight  into  the  linear  behavior  of  each  instability  in  the 
important  "large-scale''  regime. 

The  purpose  of  this  paper  is  to  derive  a  single,  general,  ccnlocal 
dispersion  equation  which  describes  the  relevant  interchange 
instabilities:  Raylelgh-Taylor,  E  x  B  gradient  drift,  and  current 
convective.  The  important  equilibrium  quantities  considered  in  the 
derivation  are  gravity,  an  ambient  electric  field,  a  current  parallel  to 
the  background  magnetic  field,  and,  of  course,  a  density  gradient.  Simple 
analytic  expressions  of  the  growth  rate  for  each  instability  in  the  long 
wavelength  limit  are  presented  in  both  the  collisional  (v^  »  u)  and 
collisionless  (u  »  viQ)  regimes,  where  u  is  the  wave  frequency  and  is 
the  ion-neutral  collision  frequency. 
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The  organization  of  the  paper  Is  as  follovs.  In  the  next  section  we 
derive  the  general  nonlocal  aode  equation.  In  Section  III  we  derive  the 
nonlocal  dispersion  equation.  In  Section  IV  we  present  simple  analytical 
expressions  for  the  growth  rates  of  the  various  interchange  instabilities 
in  the  long  wavelength  limit.  In  the  final  section  we  present  our 
conclusions . 

II.  DERIVATION  OF  MODE  EQUATION 

The  plasma  configuration  and  slab  geometry  used  In  the  analysis  are 
shown  in  Fig.  la.  The  ambient  magnetic  field  is  in  the  z  direction 

^  a 

(3  *  Bq  ez),  the  ambient  electric  field  is  in  the  y  direction  (E  »  EQ  ey), 

the  ambient  current  is  in  the  z  direction  (J  *  Jq  e^),  gravity  Is  in  the  x 

direction  (^  *  g  exj  and  the  density  is  taken  to  be  Inhomogeneous  in  the  x 

direction  (n  »  nQ(x)].  A  collislonal  plasma  is  assumed  such  that 

v  .,/fl.  «  1»  v  /G  «  1*  v.  /Q.  «  1,  and  v.  „/G.  «  1  where  3  »  eBn/m  c 
ei  e  en  e  le  1  in  1  a  0  a 

is  the  cyclotron  frequency  of  species  a,  \>ei  refers  to  electron-ion 

collisions,  u.  to  electron-neutral  collisions,  v.  to  ion-electron 
en  le 

collisions,  and  to  ion-neutral  collisions.  Furthermore,  we  assume  that 

ven^e  ^  vin^i  our  ana^78^s*  Perturbed  quantities  are  assumed  to  vary 

as  5p  -  6p(x)  exp  [i(k^y  +  k^z  -  ut)]  and  it  is  assumed  that  «  1, 

kp,  «  1,  and  k  «  k  ,  where  p.  is  the  mean  ion  Larmor  radius.  That  is, 
1  z  y  1 

we  consider  low  frequency,  long  wavelength,  field-aligned  perturbations. 
We  neglect  temperature  effects.  Finally,  we  consider  only  electrostatic 
oscillations  and  assume  quasi-neutrality  (n#  »  ^j. 


The  fundamental  aquations  used  in  the  analysis  are  continuity, 
•omentum  transfer,  and  charge  conservation.  In  the  neutral  frame  of 


reference: 


^2+  V.(nVj  -  0  (1) 
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(2) 

(3) 


7*l  m  *  0 


(*) 


where  a  denotes  species  (e:  electrons,  1:  ions)  and  ocher  variables  have 
their  usual  meaning.  Note  that  we  have  neglected  electron  inertia  effects 
in  Eq.  (2)  but  have  retained  ion  inertia  effects  in  Eq.  (3).  However,  we 
have  assumed  3/3t  »  7  «7  in  Eq.  (3)  which  is  justified  for  the 
instabilities  under  consideration#  The  equilibrium  drifts  are  given  by 


‘  0 


7  ♦  r- A.  + 

-10  %  %  ni  Bo 


e  +  7 ,  e 
y  d  z 


(5) 

(6) 


where  we  have  taken  7X  -  7X  -  cEq/B  and  have  chosen  the  ions  to  carry  the 

A 

parallel  current  (J  •  enV^  e^). 
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We  now  linearize  Eqs.  (l)-(3)  and  Cake  n  -  n«  +  Sn,  V  ■  V  n  +  5V  , 

u  -ct  s*au 

and  E  ■  -  7$  where  $  is  Che  perturbed  electrostatic  potential.  Using  Eqs. 
(2)  and  (3)  we  find  that 


‘’.■•tV  •*  ‘z 


which  can  be  written  as 
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(8) 
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where  u  »  u  +  iviQ  and  the  prime  indicates  a  derivative  with  respect  to 
x.  We  have  neglected  collisional  effects  on  the  electron  notion 
perpendicular  to  B  and  on  the  ion  notion  parallel  to  3.  This  is  Justified 
since  we  have  taken  v  / J2  <<  u/n,  which  is  appropriate  for  ionospheric 
conditions  in  the  F  region. 

We  now  substitute  Eqs.  (9)  and  (10)  into  Eq.  (4)  and  obtain 


(U) 


5 


•T^dsTa 


We  relate  $n  and  $  using  the  electron  continuity  equation  and  find  that 
(from  Eqs.  (1)  and  (9)) 


60  "  f  ^  [ao*'  '  (V)-  +  ino  r-^T"  k z*] 

y  ei 


(12) 


Finally,  substituting  Eq.  (12)  into  Eq.  (11)  we  arrive  at  the  mode  equation 
which  describes  the  interchange  instabilities  relevant  to  ionospheric 
structuring 


i  k  V.  n, 

(V')'  “  aoky ^ 1  +  1  T  v — t0« 

k  ei  aj 
7 

Q.  k  v,  cE_ 
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(13) 


The  first  term  in  Eq.  (13)  relates  to  the  node  structure  in  the  x 
direction,  the  second  term  relates  to  the  mode  structure  in  the  y  and  z 
directions,  and  the  final  term  contains  the  driving  mechanisms  for  the 
various  instabilities  (viz.,  the  Rayleigh-Taylor  instability  (*  gn'],  the 
E  x  B  gradient  drift  instability  («  EqUq),  and  the  current  convective 
instability  «  V^n^)). 


III.  ANALYSIS  OF  MODE  EQUATION 

As  noted  in  the  introduction,  the  bulk  of  linear  analyses  of  the 
Interchange  instabilities  under  study  have  made  use  of  the  local 
approximation  so  that  only  short  wavelength  modes  are  considered.  That  is. 
It  is  assumed  that  k^L^  »  k\^  »  l  where  L  •  (a in  n^/ 3x)  *  is  the  density 
gradient  scale  length  and  kx  -  3/3x.  With  this  assumption  one  can  neglect 
the  first  term  on  the  LHS  of  Eq.  (13)  and  obtain  a  relatively  simple 


dispersion  equation 


(14) 


where  r  -  (1  +  i(k2/k2 )  (fl^/w)  )  j^2  •  Equation  (14)  has  been 
thoroughly  analyzed  in  the  literature  for  each  of  the  Instabilities  under 
consideration  (see  Ossakow  et  al.  (1982)  and  references  therein)*  However, 
we  are  considering  the  long  wavelength  Unit  (kL  <<  l)  and  Eq.  (13)  mist  be 
solved  in  order  to  obtain  the  nonlocal  dispersion  equation. 

For  simplicity,  we  consider  a  density  profile  with  a  single 
discontinuity  at  x  ■  0  (see  Fig.  lb)  given  by 
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x  >  0 
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For  x  *  0,  n-*  »  0  and  Eq.  (13)  reduces  to 
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where 


r2  - 


,  + 1 < V. 

*  +  1  71  ~ — 
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the  solution  to  Eq.  (16)  is  taken  to  be 
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(16) 
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Since  the  modes  are  assumed  Co  be  bounded  as  x  +  •  we  note  chat 


$(x)  * 


rx 


k  Vx 

.  7 


x  >  0 

x  <  0 


(19) 


We  require  chat  Che  tangential  component  of  the  electric  field  be 
continuous  at  x  m  0  (Hasegawa,  1971)  which  means  that  $  is  continuous  at  x 
•  0.  This  is  equivalent  to  requiring  that  the  interface  velocity  and  the 
fluid  velocity  perpendicular  to  the  interface  be  equal  (Chandresekhar , 
1961),  i.e.,  is  continuous  at  the  discontinuity.  Thus,  »  *2  ia  Eq. 
(19)  so  that 


$(x)  - 
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-kyrx 


k  rx 
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x  >  0 


x  <  0 


(20) 


To  obtain  the  dispersion  equation  we  integrate  Eq.  (13)  across  the 
discontinuity  at  x  ■  0.  Thus,  we  have 


Av'i'41  -  /'t“ok7r2  ■^lkr!r'l‘rrr'  v,»>o>  - 0 

“t  "€  U>  i  1  0 


(21) 


Since  $  is  continuous  across  the  boundary  at  x  -  0,  it  is  found  chat  Eq. 
(21)  leads  to  (letting  e  •*>  0) 


(v')i 


in 
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-  k  V  )  ( n 
2  d  ‘  -  1 


(22) 


where  (1,2)  indicate  the  regions  x  >  0  (+s)  and  x  <  0  (-s),  respectively. 
Substituting  Eq.  (20)  into  Eq.  (22) 
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which  is  the  general  dispersion  equation  which  describes  long  wavelength 
interchange  modes. 

It  is  interesting  to  compare  the  short  wavelength  (i.e.,  local)  and 

long  wavelength  (i.e.,  nonlocal)  dispersion  aquations,  that  is,  Zqs.  (11) 

and  (23),  respectively.  These  equations  are,  in  fact,  quite  similar.  We 

define  an  effective  wavenumber  k  given  by  k  =  k  F  *  (k“  +  ik“'n,/mj 

y  y  m  . 

(Jle/vei)j  and  use  the  definition  of  the  density  gradient  scale  length 
[L  «  [n'/ng]  to  note  the  following.  The  transition  from  the  short 
wavelength  limit  to  the  long  wavelength  limit  can  be  made  by  substituting 
(n^  -  n^j/fn^  +■  n^)  for  [kLJ  ^  in  Eq.  (14).  By  doing  this  Eq.  (23)  is 
recovered.  A  significant  feature  of  this  substitution  is  that  the 

wavenumber  dependence  of  the  growth  rates  in  the  long  and  short  wavelength 
limits  will  be  different.  It  will  be  shown  that  the  short  wavelength 

growth  rate  is  independent  of  k  while  the  long  wavelength  growth  rate  Is 
proportional  to  k”P  where  p  <  1  and  the  actual  value  of  p  depends  upon  the 
plasma  parameters. 

Finally,  we  close  this  section  by  rewriting  Eq.  (23)  in  a  form  which 
is  in  line  with  ionospheric  structuring  terminology.  That  is,  we  define  M 
“  n^/n2  which  can  be  thought  of  as  the  ratio  of  Pedersen  conductivity 


9 


ins id*  eh*  plasma  cloud  to  that  outs  id*  Che  plasma  cloud.  With  this 
definition,  Eq.  (23)  takes  the  form 


vla  “*0 
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(24) 


IV.  ANALYTICAL  RESULTS 

In  this  section  we  present  analytical  expressions  for  the  growth  rates 
of  Che  interchange  instabilities  under  study  in  both  the  collisional 
(via  >>  <u)  and  collisionless  limits  (u»  »  Simplifying  assumptions 

will  be  made  in  order  to  distill  essential  parameter  scalings  of  the  growth 
rates. 

A.  Colli3ional  Limit 

In  this  limit  Z  *  ivln  and  r2  *  I  +  (k2/k2)  (q^/v^  j  { / v£a  3  *  The 
general  nonlocal  dispersion  equation  is  then  given  by 


Y  ’ 
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fli  . 
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r£ _ iE 
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(25) 


where  the  subscript  ni  denotes  nonlocal.  Instability  (i.e.,  y  >  0)  can 
occur  for  M  >  1  when  g  <  0  (Rayleigh-Taylor) ,  Eq  >  0  (E  x  B  gradient 
drift),  and  V^  >  0  (current  convective).  We  note  that  the  local  (short 
wavelength)  growth  rate  is  given  by 


Y 


i 


t\uL 


(26) 


which  is  independent  of  ky.  The  nonlocal  (long  wavelength)  growth  rate  can 
then  be  expressed  as 


Qi 


■  v  (Kr)r  ->1 


(27) 


Thus,  che  growth  rates  of  the  collisioaal  interchange  instabilities 
relevant  to  ionospheric  structuring  in  the  long  wavelength  limit  are 
directly  proportional  to  ky(5f-l)/(ltf-l).  This  limit  applies  to  the  low 
altitude  F  region  (<  400  km). 

B.  Collisionless  Limit  (w  »  ) 

In  this  limit  u>  *  u>  so  that  *  1  +  i(k^/k^)  )  (yu)  •  Since 

can  be  a  function  of  m  in  this  case  (i.e.,  for  k2  finite),  che  relationship 

between  the  local  and  nonlocal  growth  rates  is  not,  in  general,  as  simple 

as  in  the  collisional  case.  We  consider  then  two  limits:  kz  *  0  and 

fc,  *  0.  For  the  case  k_  *  0  the  relevant  instabilities  are  the  Rayleigh- 
z  z 

Taylor  and  che  E  x  3  gradient  drift.  For  kz  *  0  the  relevant  instability 
is  the  current  convective.  These  results  apply  to  the  high  altitude  F 
region  (>  500  km). 

1.  -  0 

In  the  collisionless  flute  mode  limit  (k«B  *0),  we  find 

v  w 

that  r  *  1  so  that  the  nonlocal,  long  wavelength  growth  rate  is  simply 
(from  Eq.  (24)) 


n! 


fi  _  *0,1/2,.  L  -  1 1 -1 1/ 2 


(29) 


Instability  can  occur  for  M  >  1  when  g/L  <  0  (Rayleigh-Taylor)  and  E^/L  >  0 
(E  x  B  gradient  drift).  Again,  we  note  that  the  local,  short  wavelength 
growth  rate  is  given  by 


11 


(30) 


Tt  '•L  L  Bo 


so  that 


(31) 


Thus,  in  contrast  to  the  collislonal  Rayleigh-Taylor  and  E  *  B  gradient 
drift  Instabilities,  the  growth  rate  of  the  colllslonless  long  wavelength 
mode  is  proportional  to  [kyL(>t-l)/ (»+■!)]  ^^. 

2.  s  0 

In  the  colllslonless  limit  with  k,  finite  ve  consider  only  the 

Z 

currenc  convective  Instability  (i.e.,  we  take  g  ■  0  and  Eg  *  0).  In  this 
2  2  2 

limit  r  ■  1  +  ifk  /k  ':fa  /\>  .  ]fa,  / u]  so  that  the  long  wavelength,  nonlocal 
v  z  y '  ^  e  ei J -  l 

dispersion  equation  is  given  by 


kz  ni  %  ,1/2 


«  (1  +  i  ) 

v  ,2  to  v,  ; 


ei 


M  -  1 

kz7d^i  M  +  1 


(32) 


We  obtain  a  simple  expression  for  the  growth  rate  by  assuming  that  m/a^  « 
(k^/ky) (0e/vei ).  In  this  limit  Eq.  (32)  becomes 


3 

to 


i_  a  ,_2^„  el  ,M  ■*  1,1 
+  1  yd°i  Q  +  iJ 


(33) 


We  find  that 


a.i 


m  rv*i  T,l/3,  T,2/3,M  -  1,2/3,’d, 

\im(—  T~]  (k7L)  (mTT)  C“) 


2/3, M  -  1,2/3, 


(34) 
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We  can  relate  this  nonlocal  growth  rate  to  the  maximum  local  growth 


rate  which  Is  given  by  (Chaturvedi  and  Ossakow,  1981) 


Y 


oi 


A 

4v  J 


1/3 


(35) 


Thus,  we  find  that 


rn£  "  °fV(i^T^2/3Y«i  (36> 

where  a  is  a  factor  of  order  unity.  Thus,  the  growth  rate  of  the  long 
wavelength,  collisionless  current  convective  instability  is  proportional  to 

Pc7L(M-1)/0*H)]2/3. 

V.  CONCLUSION 

We  have  derived  a  general,  nonlocal  dispersion  equation  (Eq.  (13)) 
which  deacrlbes  the  long  wavelength  Halt  (kL  «  1)  of  several  Interchange 
instabilities  (i.e.,  layleigh-Taylor,  E  x  B  gradient  drift,  and  current 
convective).  Analytical  expressions  for  the  growth  rates  of  these 
instabilities  have  been  presented  in  both  the  collislonal  (\>ia  »  w)  and 
collisionless  (a  »  )  regimes.  An  1 sport ant  result  is  the  scaling 

of  with  kyL  and  ti.  These  scalings  are  summarized  in  Table  I.  The 
implications  of  this  work  in  regard  to  ionospheric  structuring  are 
discussed  by  Zalesak  and  Hubs,  and  by  Wlttwer  in  these  proceedings. 
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Table  1  —  Scaling  of  growth  rates  with  respect  to  k  and  M 

Instability  Colllslonal  Collisionless 

(uia  >>  *0  >>  vin) 


Rayleigh-Taylor 

V$T 

4) 

tkyL^M  + 

i)]1/2 

E  x  B  gradient  drift 

V  v  w 

vfir 

t) 

[kyL(}Tir 

Current  convective 

k  L(-~- 
y  +■ 

*> 

t)’2/3 

IS 
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HUES 

COURT  STREET 

SYRACUSE,  ST  13201  , 

OICT  AITS  a.  MILLNAN 

GEOPHYSICAL  XSSTtTDTS 

university  or  alaska 
ruums,  ax  99701 

(ALL  CLASS  ATTN:  SECT* ITT  OFFICEE) 
OICT  ATTN  T.S.  DAVIS  (UNCLASS  ONLY) 
OICT  ATTN  TECHNICAL  LIBRARY. 

OICT  ATTN  SEAL  BROWN  (UNCLASS  ONLY) 

GTE  STL  TANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTKRH  01 7 
77  A  STKEST 
NEEDHAM,  MA  0219* 

OICT  ATTN  DICK  STEINBOP 

ass,  inc. 

2  altoed  circle 

SEDFOtD,  MA  01730 

OICT  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAICN,  CL  41820 

(ALL  CORSES  ATTN  DAN  MCCLELLAND) 

OICT  ATTN  L  TEH 

INST  iron  FOR  DEFENSE  ANALTSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA.  VA  22311 
OICT  ATTN  J.H.  AEZN 
OICT  ATTN  ERNEST  BADER 
OICT  ATTN  BANS  WOLFAIS 
OICT  ATTN  JOEL  BENGSTON 

INTL  TEL  A  TELEGRAPH  CORPORATION 
SOO  WASHINGTON  AVENUE 
NUTLET,  NJ  07110 

OICT  ATTN  TECHNICAL  LIBRARY 

JATCOR 

11011  TOIREYANA  BOAS 
F.O.  BOX  8313* 

SAM  DIEGO,  CA  92138 

OICT  ATTN  J.L.  SPERLING 


JOHNS  BOPEZNS  UNIVERSITY 
APPLIED  PHYSICS  LARORATORT 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OICT  ATTN  DOCUMENT  LIBRARIAN 
OICT  ATTN  THOMAS  PQTEMRA 
OICT  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORF 
P.O.  30X  74*3 

COLORADO  SPRINGS,  CO  *0933 
OICT  ATTN  T.  MEAGHER 

RAMAN  TEMPO-COTTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICT  ATTN  DAS  LAC 
OICT  ATTN  WARREN  S.  KNAPP 
OICT  ATTN  WILLIAM  MCNAMARA 
OICT  ATTN  B.  OAKRILL 

LINEAR  IT  CORP 
10453  ROSELLE 
SAN  DI2G0,  CA  92121 
OICT  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  ISC 
P.O.  SOX  504 
SUNNYVALE,  CA  94088 
OICT  ATTN  DEPT  60-12 
OICT  ATTN  D.2.  CHURCHILL 

LOCKHEED  MISSILES  *  SPACE  CO.,  INC. 

3231  HANOVER  STREET 
PALO  ALTO,  CA  9A304 

OICT  ATTN  MARTIN  WALT  DEPT  52-12 
OICT  ATTN  W.L.  IMHOF  DEFT  52-12 
OICT  ATTN  RICHARD  C.  JOHNSON  DEPT  52-12 
OICT  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32305 
OICT  ATTN  R.  HEFFNIR 

M.I.T.  LINCOLN  LABORATORY 
F.O.  BOX  73 
LXXINGTON,  MA  02173 

OICT  ATTN  DAVID  M.  TOWLE 
OICT  ATTN  L.  L0UGHL2N 
OICT  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
3301  30LSA  AVENUE 
aUNTINCTON  REACH,  CA  92647 
01CT  ATTN  5.  HARRIS 
01CT  ATTN  J.  HOULE 
01CT  ATTN  GEORGE  MXOZ 
OICT  ATTN  '4.  OLSON 
OXCT  ATTN  a.W.  HALPRIN 
OICT  ATTN  TECHNICAL  LIS EAST  SERVICES 

MISSION  RESEARCH  CORPORATION 
733  STATE  STREET 
SANTA  RARSARA,  CA  93101 
OICT  ATTN  P.  PISCHER 
OICT  ATTN  U.P.  CREVIER 
OICT  ATTN  STEVEN  L.  GUTSCHZ 
OICT  ATTN  R.  BOGUSCH 
OICT  ATTN  R.  HENDRICK 
OICT  AITS  RALPH  SIL3 
OICT  ATTN  DAVE  SOWLE 
OICT  ATTN  ?.  FAJEN 
OICT  ATTN  M.  SCHEIBE 
OICT  ACTS  CONRAD  L.  LONGMIRZ 
OICT  ACTS  S.  WHITE 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.Z. 

ALBUQUERQUE,  NEW  MEXICO  87106 
OICT  R.  STELLINGWERF 
OICT  M.  ALNE 
OICT  L.  BRIGHT 

MITRE  CORPORATION.  THE 
P-0*  BOX  208 
BEDFORD,  HA  01730 

OICT  ATTN  JOHN  MORGAN  STERN 
OICT  ATTN  G.  HARDING 
OICT  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

VESICATE  RESEARCH  PARK 
1820  QOLLT  MADISON  BL7D 
MCLEAN,  VX  22101 
OICT  ATTN  V.  HALL 
OICT  ATTN  V.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90023 
OICT  ACTS  E.C.  FIELD,  JR. 


PENNSTL  VANIA  STATS  UNI  VERS ITT 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENCISEERINC  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICT  ACTS  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  ISC. 

4  ARROW  DRIVE 
aOBCRN,  MA  01801 

OICT  ACTS  IRVING  L.  CDF  SKI 

PHYSICAL  DYNAMICS,  ISC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICT  ACTS  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  RET,  CA  90291 
OICT  ATTN  FORREST  GILMORE 
OICT  ATTN  WILLIAM  3.  WRIGHT,  JR. 
OICT  ATTN  ROBERT  ?.  LZLE  VIES 
OICT  ATTN  WILLIAM  J.  RARZAS 
OICT  ATTN  H.  ORT 
OICT  ATTN  C.  MACDONALD 
OICT  ATTN  R.  TURCO 
OICT  ATTN  L.  Da RAND 
OICT  ATTN  W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICT  ATTN  CULLEN  CRAIN 
OICT  ATTN  ED  3EDR0ZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 
OICT  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  VEST  42nd  STREET 
NEW  YORK,  ST  10036 

OICT  ACTS  VINCE  TRAPANI 


1 


SCIENCE 

APPLICATIONS ,  ESC. 

TRW  DEFENSE  &  SPACE  STS  CROC? 

1130  PROSPECT  PLAZA 

ONE  SPACE  PARE 

LA  JOLLA.  CA  92037 

REDONDO 

BEACH,  a  90278 

OICT 

ATTH  LEWIS  X.  LIHSOH 

OICT 

ATTH  R.  E.  PLEBUCZ 

OICT 

ATTH  DANIEL  A.  HAMLIN 

OICT 

ATTH  S.  ALTSCHULER 

OICT 

ATTH  E.  FREEMAN 

OICT 

ATTH  D.  DEE 

OICT 

ATTH  S.A.  STIAEER 

OICT 

ATTH  D/  STOC  EWELL 

OICT 

ATTH  CURTIS  A.  SMITH 

SNTT/1573 

OICT 

ATTH  JACK  MCDOOGALL 

7ISIDTNE 

SCXESCS 

APPLICATIONS ,  ESC 

SOUTH  BEDFORD  STREET 

1710  COODRIDGR  A. 

BUILXSCTON,  HASS  01803 

MCLEAN, 

VA.  22102 

OICT 

ATTH  W.  REIDT 

ATTH: 

J.  C0C1CATHE 

OICT 

ATTH  J.  CARPEHTER 

OICT 

ATTH  C.  HUMPHRET 

SU  INTERNATIONAL 
333  RAVEHSUOOD  AVEHUE 
MENLO  PARK,  CA  9*023 


OICT 

ATTH 

DONALD  HKXLSON 

OICT 

ATTH 

ALAN  BURNS 

OICT 

ATTH 

0.  SMITH 

OICT 

ATTH 

R.  TSUNOOA 

OICT 

ATTH 

DAVID  A.  JOHNSON 

OICT 

ATTH 

Walter  >.*  cmsHUT 

OICT 

ATTH 

CHARLES  L.  UNO 

OICT 

ATTH 

WALTER  JATE 

OICT 

ATTH 

J.  VICKRET 

OICT 

ATTH 

RAT  L.  LZADABRAND 

OICT 

ATTH 

S.  CARPENTER 

OICT 

ATTH 

C.  PRICE 

OICT 

ATTH 

R.  LIVINGSTON 

OICT 

ATTH 

V.  GONZALES 

OICT 

ATTH 

D.  MCDANIEL 

techholoct  nrrEuunouL  an? 
73  WIGGINS  AVENUE 

Bedford,  ma  01730 

OICT  ATTN  W.P.  BOqUIST 


CO. 


TOTON  RESEARCH 
P.0.  Sox  6890 
SANTA  BARBARA,  CA  93111 

oicr  an  johm  is*,  a. 

OICT  ATTH  JOEL  GARBARIHO 
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END 

'DATE 

FILMED 


DTIC 


